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This paper addresses the symbolic representation of non-convex real polyhedra, i.e., sets of real vec- 
tors satisfying arbitrary Boolean combinations of linear constraints. We develop an original data 
structure for representing such sets, based on an implicit and concise encoding of a known structure, 
the Real Vector Automaton. The resulting formalism provides a canonical representation of poly- 
hedra, is closed under Boolean operators, and admits an efficient decision procedure for testing the 
membership of a vector. 



1 Introduction 

Algorithms and data structures for handling systems of linear constraints are extensively used in many 
areas of computer science such as computational geometry lH3l . optimization theory ll24ll . computer- 
aided verification ||T0llT5l , and constraint programming 11231 . In this paper, we consider systems defined 
by arbitrary finite Boolean combinations of linear constraints over real vectors. Intuitively, a non-trivial 
linear constraint in the ^-dimensional space describes either a (n — 1) -plane, or a half-space bounded by 
such a plane. A Boolean combination of constraints thus defines a region of space delimited by planar 
boundaries, that is, a polyhedron (also called n-polytope). 

Our goal is to develop an efficient data structure for representing arbitrary polyhedra, as well as asso- 
ciated manipulation algorithms. Among the requirements, one should be able to build representations of 
elementary polyhedra (such as the set of solutions of individual constraints), to apply Boolean operators 
in order to combine polyhedra, and to test their equality, inclusion, emptiness, and whether a given point 
belongs or not to a polyhedron. 

A typical application consists in representing objects in a 3D modeling tool, in which shapes are 
approximated by polyhedral meshes. By applying Boolean operators, the user can modify an object, for 
instance, drilling a circular hole amounts to computing the Boolean difference between the object and a 
polyhedron approximating a cylinder. This application requires an efficient implementation of Boolean 
operations: A local modification performed on a complex object should ideally only affect a small part 
of its representation. 

Another application (actually our primary motivation for studying this problem) is the symbolic 
representation of the reachable data values computed during the state-space exploration of programs. 
In this setting, a reachable set is computed iteratively, by repeatedly adding new sets of values to an 
initial set, and termination is detected by checking that the result of an exploration step is included 
in the set of values that have already been obtained. In this application, it is highly desirable for a 
representation of a set to be independent from the history of its construction, since reachable sets often 
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have simple structures, but are computed as the result of long sequences of operations. We are particularly 
interested in linear hybrid systems Q, for which symbolic state-space exploration algorithms have been 
developed CD COD, requiring efficient data structures for representing and manipulating systems of linear 
constraints. Existing representations either fail to be canonical [ 15 , 17 ], or impose undue restrictions on 
the linear constraints that can be handled ifTTTl . 

For some restricted classes of systems of linear constraints, data structures with good properties are 
already well known. Consider for instance conjunctions of linear constraints, which correspond to convex 
polyhedra. A convex polyhedron can indifferently be represented by a list of its bounding constraints, or 
by a finite set of vectors (its so-called vertices and extremal rays) that precisely characterize its shape (221. 
An efficiently manageable representation is obtained by combining the bounding constraints and the 
vertices and rays of a polyhedron into a single structure |[T0l[T9l l3l. 

There are several ways of obtaining a representation suited for arbitrary combinations of linear con- 
straints. A first one is to represent a set by a logical formula in additive real arithmetic. This approach 
is not efficient enough for our intended applications, since testing set emptiness, equality, or inclusion 
become NP-hard problems [12J. A second strategy is to decompose a non-convex polyhedron into an 
explicit union of convex polyhedra (which may optionally be required to be pairwise disjoint). The main 
disadvantage of this method is that a set can generally be decomposed in several different ways, and that 
checking whether two decompositions correspond to the same set is costly. Moreover, simplifying a long 
list of convex polyhedra into an equivalent shorter union is a difficult operation. 

Another solution is to use automata [8]. The idea is to encode ^-dimensional vectors as words over 
a given alphabet, and to represent a set of vectors by a finite-state machine that accepts the language of 
their encodings. This technique presents several advantages. First, with some precautions, computing 
Boolean combinations of sets reduces to applying the same operators to the languages accepted by the 
automata that represent them, which is algorithmic ally simple. Second, provided that one employs de- 
terministic automata, checking whether a given vector belongs to a set becomes very efficient, since it 
amounts to following a single path in a transition graph. Finally, some classes of automata can easily be 
minimized into a canonical form. This approach has already been applied successfully to the represen- 
tation of arbitrary combinations of linear constraints, yielding a data structure known as the Real Vector 
Automaton (RVA) JUIU. 

Even though RVA provide a canonical representation of polyhedra, and admit efficient algorithms for 
applying Boolean operators, they also have major drawbacks. First, they cannot handle efficiently linear 
constraints with coefficients that are not restricted to small values, since the size of RVA generally gets 
proportional to the product of the absolute values of these coefficients [9]. Second, RVA representing 
subsets of the 7i-dimensional space get unnecessarily large for large values of n. 

The contribution of this paper is to tackle the first drawback. We introduce a data structure, the 
Implicit Real Vector Automaton (IRVA), that represents polyhedra in a functionally similar way to RVA, 
but much more concisely. The idea is to identify in the transition relation of RVA structures that can be 
described efficiently in algebraic notation, and to replace these structures by their implicit representation. 
We show that checking whether a vector belongs to a set represented by an IRVA can be decided very 
efficiently, by following a single path in its transition graph. We also develop algorithms for minimizing 
an IRVA into a canonical form, and for applying Boolean operators to IRVA. 
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2 Basic Notions 

2.1 Linear Constraints and Polyhedra 

Let n G N>o be a dimension. A linear constraint over vectors x G W is a constraint of the form a.x#b, 
with a G Z", fceZ, and # G {<,<,=, >, >}. A finite Boolean combination of such constraints forms 
a polyhedron. If a polyhedron can be expressed as a finite conjunction of linear constraints, it is said to 
be convex. A polyhedron that can be expressed as a conjunction of linear equalities, i.e., constraints of 
the form a.x = b, is an affine space. An affine space that contains is a vector space. The dimension 
dim(VS) of a vector space VS is the size of the largest set of linearly independent vectors it contains. 

Finally, given a convex polyhedron D C W, a polyhedron PCS", and a vector v£D, we say that P 
is conical in D with respect to the apex v iff for all x G D and A G]0, 1 [, we have * G P 44> A (x — v) + v G P. 
(Intuitively, this condition expresses that within D, the polyhedron P is not affected by a scaling centered 
on v.) It is shown in [7 ] that the set of the vectors v with respect to which P is conical in D necessarily 
coincides with an affine space over D. 

2.2 Real Vector Automata 

This section is adapted from (HHHTl. Let r G N>i be a numeration base. In the positional number system 
in base r, a number z G M>o can be encoded by an infinite word a p -\a p -2 ■ ■ ■ ao *a_ia_2fl'-3 • • ., where 
V/ : a, G {0, 1, . . . ,r — 1}, such that z = T,i<p a i r '- (The distinguished symbol "★" separates the integer 
from the fractional part of the encoding.) Negative numbers are encoded by using the r's-complement 
method, which amounts to representing a number z G R<o by the encoding of z + r p , where p is the length 
of its integer part. This length p does not have to be fixed, but must be large enough for the constraint 
— r p ~ x < z < r p ~ l to hold, in order to reliably discriminate the sign of encoded numbers. Under this 
scheme, every real number admits an infinite number of encodings in base r. Note that some numbers 
admit different encodings with the same integer-part length, for instance, the base-2 encodings of 1/4 
form the language + *010 ffl U + ★001 £0 . Such encodings are then called dual. 

The positional encoding of numbers generalizes to vectors in W 1 , with n G N>o- A vector is en- 
coded by first choosing encodings of its components that share the same integer-part length. Then, these 
component encodings are combined by repeatedly and synchronously reading one symbol in each com- 
ponent. The result takes the form of an infinite word over the alphabet {0,1,.../-1}"U {*} (since the 
separator is read simultaneously in all components, it can be denoted by a single symbol). It is also worth 
mentioning that the exponential size of the alphabet can be avoided if needed by serializing the symbols, 
i.e., reading the components of each symbol sequentially in a fixed order rather than simultaneously 0. 

This encoding scheme maps any set S C W onto a language of infinite words. If this language is 
a)-regular, then it can be accepted by an infinite-word automaton, which is then known as a Real Vector 
Automaton (RVA) representing the set S. 

Some classes of infinite- word automata are notoriously difficult to handle algorithmically ll25l . A 
weak automaton is a Biichi automaton such that each strongly connected component of its transition 
graph contains either only accepting or only non-accepting states. The advantage of this restriction is 
that weak automata admit efficient manipulation algorithms, comparable in cost to those suited for finite- 
word automata [26]. The following result is established in j8). 

Theorem 1 Let n G N>o- Every polyhedron of W can be represented by a weak deterministic RVA, in 
every base r G N>i. 
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In the sequel, we will only consider weak and deterministic RVA. These structures can efficiently be 
minimized into a canonical form [21], and combining them by Boolean operators amounts to performing 
similar operations on the languages they accept. Implementations of RVA are available as parts of the 
tools LASH ED and LIRA HI. 

3 The Structure of Polyhedra 

It is known that RVA can form unnecessarily large representations of polyhedra. For instance, a finite- 
state automaton recognizing the set of solutions {x\,X2) of the constraint x\ = r k X2 in base r essentially 
has to check that x\ and X2 have identical encodings up to a shift by k symbols, and thus needs 
states for its memory. On the other hand, the algebraic description of the constraint x\ = r k X2 requires 
only 0(k) symbols. 

In this section, we study the transition relation of RVA representing polyhedra, with the aim of finding 
internal structures that can more efficiently be described in algebraic notation. 

3.1 Conical Sets 

It has been observed in [ 14 ] that, for every polyhedron P CW and point vGl", the set P is conical in all 
sufficiently small convex neighborhoods of v. We now formalize this property, and prove it by reasoning 
about the structure of RVA representing P. This will provide valuable insight into the principles of 
operation of automata-based representations of polyhedra. 

For every v = (vi,V2,. . . ,v«) 6l" and £ G M>o, let Af e (v) denote the rc-cube of size [e]" centered on 
v, that is, the set [vi — e/2,v\ +e/2] x [v 2 — e/2,v 2 + e/2] x • • • x [v„ — e/2,v„ + e/2]. 

Theorem 2 Let P C W be a polyhedron, with n G N>o, and let v 6l" be an arbitrary point. For every 
sufficiently small e G M>o, the set P is conical in N e (v) with respect to the apex v. 

Proof: Let sz? be a RVA representing P in a base r G N>i, which exists thanks to Theorem [T] We 
assume w.l.o.g. that srf is weak, deterministic, and has a complete transition relation. Consider a word w 
encoding v in base r. For each k G N, let w\ denote the finite prefix of w with k symbols in its fractional 
part, i.e., such that Wk = u*u' with \u'\ = k. The set of all vectors that admit an encoding of prefix 
Wk forms a «-cube C Wk of size [r~ k ] n . For every k G N, we have v G C Wk and C Wk+l C C Wk , leading to 
PLteN^W = {^}- Intuitively, each symbol read by srf reduces by a factor r" the size of the set of possibly 
recognized vectors. 

Consider e G R>o with e < 1. The set A^(v) is covered by the union of the sets C Wk for all w 
encoding v, choosing k such that r~ k > e. It is thus sufficient to prove that for every word w encoding 
v and sufficiently large k, the set P is conical in C Wk with respect to the apex v. This property has been 
proved in [7 ], where it is additionally shown that the suitable values of k include those for which w^ 
reaches the last strongly connected component of &f visited by w. □ 

In the previous proof, the strongly connected components of turn out to be connected to conical 
structures present in P. This can be explained as follows. Consider two finite prefixes Wk and w^+d of 
w, with d > 0, such that w^+d only differs from w^ by additional iterations of cycles in the last strongly 
connected component of £/ visited by w. Since both Wk and Wk+d lead to the same state of £/, the sets 
of suffixes that can be appended to them so as to obtain words accepted by srf are identical. In order 
to be able to compare such sets of suffixes, we introduce the following notation. For each k G N, let 
c Wk = (c Wk ,i,c Wkj 2, ■ ■ ■ ,c WkM ) denote the vector encoded by W)t(0") ffl , in other words the vector such that 
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C Wk = [c Wjt) i,c Wjti i+r k ] x [c Wk:2 ,c Wk 2 + r k ] x ■ ■ ■ x [c Wjki „, c Wti „ + r *]. Given a «-cube C C 1" of size 
[A]" and a vector c € W, we then define the normalized view of P with respect to C and c as the set 
P[C, c] = ( 1 / X ) ( (P n C) — c) . In other words, this normalized view is obtained by a translation bringing 
c onto the origin 0, followed by a scaling that makes the size of the rc-cube in which P is observed become 
equal to [1]". 

Observe that the set P[C Wk ,c Wk ] is precisely characterized by the language accepted from the state 
of s/ reached by w&. Since this state is identical to the one reached by Wk+d, we obtain P[C Wk ,c Wk ] = 
P[Cwk+cnCwk+<i\- R eca ll that we have v E C Wk and C Wk+d C C Wk . The previous property shows that P is 
self-similar in the vicinity of v: Following additional cycles in the last strongly connected component 
visited by w amounts to increasing the "zoom level" at which the set P is viewed close to v, without 
influencing this view. It is shown in Q that this self-similarity entails the conical structure of P around 
v, which intuitively means that the zoom levels that preserve the local structure of P are not restricted to 
integer powers of r rf . 

In addition, we have established that the structure of P in a small neighborhood N £ (v) of v is uniquely 
determined by the state of srf reached by Since there are only finitely many such states, we have the 
following result. 

Theorem 3 Let P C W be a polyhedron, with n E N>o- There exists £ E M>o such that over all points 
v E W, the sets P[N e (v) , v] take a finite number of different values. Moreover, each of these sets is conical 
in [—1/2, 1 /2} n with respect to the apex 0. 

Proof: The proof follows the same lines as the one of Theorem [2] Let be a weak, deterministic, and 
complete RVA representing P in a base r E N>i. To every word w encoding a given vector v in base r, 
we associate the integer k(w) such that the path of £/ recognizing w reads the finite prefix wu w \ before 
reaching the last strongly connected component that it visits. From the previous developments, we have 
that P is conical in C Wk{w) with respect to the apex v. Furthermore, the set P [C Wjfc( . , c Wk , . ] only depends on 
the state of srf reached after reading wu w \, which are in finite number. It follows that, in arbitrarily small 
neighborhoods of v, the polyhedron P has a conical structure with respect to the apex v, and that there 
are only finitely many such structures over all vectors v. □ 



3.2 Polyhedral Components 

Theorem [3] shows that a polyhedron PCM" partitions W 1 into finitely many equivalence classes, each 
of which corresponds to a unique conical set in the «-cube [—1/2, 1 /2] n with respect to the apex 0. For 
each v E M", let C [—1/2, 1 /2] n denote the conical set associated to v by P. We call P^ the component 



of P associated to v. Recall that, as discussed in Section 2.1 the set of apexes according to which P^ is 
conical coincides with a vector space over [—1/2, 1 /2]". The dimension dim(P^) of the component Pp is 
defined as the dimension of this vector space. Finally, we say that a component P^ is in if v E P, and out 
if v P. 

An example is given in Figure [I] The triangle x\ > 1 A X2 < 2 A x\ — X2 < 1 in M 2 has three com- 
ponents of dimension corresponding to its vertices (1,0) (in), (1,2) (out) and (3,2) (out), three com- 
ponents of dimension 1 associated to its sides (two in and one out), and two components of dimension 2 
corresponding to its interior (in) and exterior (out) points. 



3.3 Incidence Relation 



In Section 3. 1 we have established a link between the components of a polyhedron PCM" and the 
strongly connected components (SCC) of a RVA representing P. We know that there exists a hierarchy 
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(a) (b) (c) 



Figure 1 : Example of (a) polyhedron, (b) polyhedral components, and (c) incidence relation. 

between the SCC of an automaton: That a SCC S2 is reachable from a SCC Si implies that every finite 
prefix that reaches a state of Si can be followed by a suffix that ends up visiting S2, while the reciprocal 
property does not hold. In a similar way, we can define an incidence relation between the components of 
a polyhedron. 

Definition 1 Let Qi , £?2 be distinct components of a polyhedron P C MP, with n G N>o- The component 
Q2 is incident to Q\, denoted Qi -< Q2, iff far all vi G W such that P$ = Qi and £ G M>o. there exists 
V2 G W such that P^ 2 = Q2 and \vi — v%\ < £■ 

Remark that the incidence relation between the components of a polyhedron is a partial order, and 
that <2i -< Qi implies dim(<2i) < dim(<22). As an example, in the triangle depicted in Figure [TJ each 
side is incident to the vertices it links, since every neighborhood of a vertex contains points from its 
adjacent sides. The reverse property does not hold. The interior and exterior components of the triangle 
are incident to each of its sides and vertices. 

3.4 How RVA Recognize Vectors 

We are now able to explain the mechanism employed by a RVA srf in order to check whether the vector 
encoded by a word w belongs or not to a polyhedron PCM". After reading an integer part and a 
separator symbol, the word w follows some transitions in the fractional part of si ', reaching a first non- 
trivial strongly connected component Si (that is, a component containing at least one cycle). At this 
location in w, inserting arbitrary iterations of cycles within Si would not affect the accepting status of w. 
This intuitively means that the prefix wt of w read so far has led us to a point that belongs to a component 
Qi of P, and that the decision can now be carried out further in an arbitrarily small neighborhood of 
this point. Reading additional symbols from w, one either stays within Si, or follows transitions that 
eventually lead to another non-trivial strongly connected component S2. Once again, this means that the 
decision can now take place in an arbitrarily small neighborhood of a point belonging to a component Q2 
of P, such that either Qi = Q2 or Qi -< Q2. The same procedure repeats itself until w reaches a strongly 
connected component that it does not leave anymore. 

In other words, in order to decide whether to accept or not a word w, the RVA si first chooses 
deterministically a component <2i of P in the vicinity of which this decision can be carried out. Then, it 
checks whether the vector v encoded by w belongs or not to Q\. If yes, the decision is taken according to 
whether <2i is in or out. If no, the RVA chooses deterministically a component Q2 incident to Q\, from 
which the same procedure is then repeated. 
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Let us now study more finely the mechanism used for moving from a component £?i that does not 
contain the vector v to another component Q2 from which v6P can be decided. One follows a path of 
s4 that leaves a strongly connected component associated to Q\, travels through an acyclic structure of 
transitions, and finally reaches a SCC associated to Q2. Recall that, as discussed in Section |3TT| at each 
step in this path, the prefix wt of w read so far determines a rc-cube C Wk . This «-cube covers some subset 
U Wk = {Pu I u G C Wk } of the components of P. If U Wk contains a single minimal component with respect 
to the incidence order -<, then this component is necessarily equal to Q2, and its associated SCC is the 
only possible destination of Wk- Indeed, all components in U Wk are then either equal or incident to Q2. If, 
on the other hand, U Wk contains more than one minimal component, then further transitions have to be 
followed in order to discriminate between them. 



4 Implicit Real Vector Automata 

Our goal is to define a data structure representing a polyhedron PCI" that is more concise than a 
RVA, but from which one can decide v£P using a similar procedure to the one outlined in Section [374] 
There are essentially three operations to consider: Selecting from a vector v an initial polyhedral com- 
ponent from which the decision can be started, checking whether v belongs or not to a given component, 
and moving from a component that does not contain v to another one from which the decision can be 
continued. We study separately each of these problems in the three following sections. 

4.1 Choosing an Initial Component 

An easy way of managing the choice of an initial component is to consider only polyhedra in which this 
component is unique. This can be done without loss of generality thanks to the following definition. 

Definition 2 Let P CM" be a polyhedron, with n G N>o- The representing cone of P is the polyhedron 
P C R" +1 = {A(xi,. . . ,x„, 1) I A G M>o A (xi, . . . ,x n ) g P}. 

For every polyhedron PCS", the polyhedron P is conical in M" +1 with respect to the apex 0, from 
which it can be inferred that every neighborhood of contains a unique minimal component £2o with 
respect to the incidence order -<. It follows that for every v G , the decision vGP can be started from 
<2o- Remark that P describes P without ambiguity, since P can be reconstructed from P by computing its 
intersection with the constraint x n+ \ = 1, and projecting the result over the first n vector components. In 
the sequel, we assume w.l.o.g. that the polyhedra that we consider are conical with respect to the apex 0. 
A similar mechanism is employed in |fl9l . 



4.2 Deciding Membership in a Component 



3.2 



Consider a polyhedron PCM" that is conical with respect to the apex 0. As explained in Section 
component of such a polyhedron is characterized by a vector space, a Boolean polarity (either in or out), 
and its incident components. Checking whether a given vector v G W belongs or not to the component 
reduces to deciding whether v belongs to its associated vector space. This is a simple algebraic operation 
if, for instance, the vector space is represented by a vector basis {b\,b2,---,b m }: One simply has to 
check whether v is linearly dependent with \b\,bz,. . . ,b m }. This approach leads to a much more concise 
representation of polyhedral components than the one used in RVA. 
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4.3 Moving from a Component to Another 

We now address the problem of leaving a component <2i of a polyhedron P C W that does not contain 
a vector v € M", and moving to a component £22 that is incident to Qi, and from which v G P can be 
decided. 

A first solution would be to borrow from a RVA representing P the acyclic structure of transitions 
leaving the strongly connected components Si associated to Q\. However, this would negate the advan- 
tage in conciseness obtained in Section |4~2) since this acyclic structure of transitions is generally as large 
as Si itself. 

The solution we propose consists in performing a variable change operation. Let {yi,yi, ■ ■ ■ ,y*m}, 
with < m < n, be a basis of the vector space associated with the component <2i- If m = n, then Q\ is 
universal and there is no possibility of leaving it. If m < n, then we introduce n — m additional vectors 
Z\, Zi, ■ ■ ■ , Zn-m, such that {y\, . . . ,y m ,zi, . . . ,z n -m} forms a basis of R". These additional vectors can be 
chosen in a canonical way by selecting among ( 1 , 0, . . . , 0) , (0, 1 , . . . , 0) , . . . , (0, 0, . . . , 1), considered in 
that order, n — m vectors that are linearly independent with {yi,y>2, ■ ■ ■ ,y m }. 

We then express the vector v in the coordinate system {yi, . . . ,y m ,Zi, ■ ■ ■ , z n -m}, obtaining a vector 
(yi,... ,y m ,Zi, ■ ■ ■ ,z n -m)- That v leaves Q\ simply means that we have (zi,... ,z n -m) / 0. As a conse- 
quence, we associate Q\ with an acyclic structure 3l\ of outgoing transitions, recognizing prefixes of 
encodings of non-zero vectors (zi, . . . ,z n -m), in order to map these vectors to the polyhedral components 
(incident to £2 1) to which they lead. 

A difficulty is that, from Theorem|2j the set P has a conical structure in arbitrary small neighborhoods 
of points in £2i- If follows that the structure S>\ has to map onto the same polyhedral component two 
vectors z and z* such that z 1 = Xz for some X G M>o- An efficient solution is to normalize the vectors 
handled by 3>\\ Given a vector z = (zi,...,z n -m) sucn that z / 0, we define its normalized form as 
[z] = (1/ (2.max;|zi|))z. In other words, [z\ is obtained by turning z into the half-line {Xz \ X G M>o}, and 
computing the intersection of this half-line with the faces of the normalization cube [—1/2,1 /2] n ~ m . In 
this way, two vectors that only differ by a positive factor share the same normalized form, and will thus 
be handled identically. 

The purpose of the structure &\ is thus to recognize normalized forms of vectors, and map them 
onto the polyhedral components to which they lead. In order to define the transition graph of Q\, 
one therefore needs a suitable encoding for normalized forms of vectors. Using the standard posi- 
tional encoding of vectors in a base r G N>i is possible, but inefficient. We instead use the follow- 
ing scheme. An encoding of a normalized vector [v\ = ([v]i, [v\2 ... , [v]„_ m ) starts with a leading sym- 
bol a G {— 2,+2, — (n — m), +(n — m)} that identifies the face of the normalization cube 
[-1/2, l/2]"~ m to which [v] belongs: If a = -i, with 1 < i <n-m, then [v],- = -1/2; if a = +i, 
then [v], = +1/2. This prefix is followed by a suffix w G {0, l} ffl that encodes the position of [v] 
within the face of the normalization cube defined by a. This suffix is obtained as follows. Assume 
that we have a G {—/,+/}, with 1 < i < n — m (which implies [v], G { — 1/2,1/2}). We turn [v\ into 
[[v]] = ([vi], . . . , [v,--i], [v,-+i], . . . , [v„_ m ]) + (1/2, 1/2, . . . ,1/2), i.e., we remove the j-th vector compo- 
nent, and offset the result in order to obtain [[v\] G [0, We then define w G {0, l} ffl as a word 
such that 0-kw is a serialized binary encoding of [[v\]. Note that some vectors v may belong to several 
faces of the normalization cube, hence their normalized form may admit multiple encodings. This is not 
problematic, provided that the structure &\ handles these encodings consistently. 

In summary, the structure 3>i is an acyclic decision graph that partitions the space of normalized 
vectors according to their destination components. Each prefix of length k read by @\ corresponds to 
a convex region R Wk C W that is conical in every neighborhood of any element of Q\, with this element 
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as apex. The situation is similar to that discussed in Section 3.4 If in a sufficiently small neighborhood 
of any point of Q\, the set of components of P covered by R Wk contains a unique minimal component Q2 
with respect to the incidence order -<, then leads to Q2. Otherwise, the decision process is not yet 
complete, and additional transitions have to be followed in !2)\ . 

4.4 Data Structure 

We are now ready to describe our proposed data structure for representing arbitrary polyhedra of W, 
with n G N>o- Recall that we assume w.l.o.g. that the polyhedra we consider are conical in W with 
respect to the apex 0. 

4.4.1 Syntax 

Definition 3 An Implicit Real Vector Automaton (IRVA) is a tuple (h,Si,Se, so, A), where 

• n is a dimension. 

• Si is a set of implicit states. Each s £ Si is associated with a vector space VS(s) C W\ and a 
Boolean polarity pol(s) G {in, out}. 

• Se is a set of explicit states, such that Se H Sj = 0. 

• so G Si is the initial state. 

• A : 5/ x ±N>o U Se x {0, 1} — > (Si U Se) is a (partial) transition relation. 

In order to be well formed, an IRVA (u,Si,Se,so,A) representing a polyhedron PCI" has to satisfy 
some integrity constraints. In particular, the transition relation A must be acyclic, and for all S\,S2 G Si 
such that A directly or transitively leads from s\ to S2, one must have VS(si) C VS^). The transition 
relation A is required to be complete, in the sense that, for every implicit state s G Si and i G N>o, A(s, +i) 
and A(s, —i) are defined iff i <n — dim(VS(s)). Furthermore, for every explicit state s G Se, both A(s,0) 
and A(s, 1) must be defined. Finally, each component of P must be described by a state in Si, and for 
every pair Q\ , Q2 of components of P such that Q\ -< Q2, there must exist a sequence of transitions in A 
leading from the implicit state associated to Q\ to the one associated to Q2. In other words, the order -< 
between the components of P can straightforwardly be recovered from the reachability relation between 
the implicit states representing them. 



4.4.2 Semantics 

The semantics of IRVA is defined by the following procedure, that decides whether a given vector v G W 
belongs or not to the polyhedron P represented by an IRVA (n,Si,SE,so,A). The principles of this 
procedure have already been outlined in Sections |4~2| and |4~3] 

One starts at the implicit state sq. At each visited implicit state s, one first decides whether v G VS(s). 
In case of a positive answer, the procedure concludes that v G P if pol(s) = in, and that v^P otherwise. In 
the negative case, the decision has to be carried out further. The vector v is transformed into v 1 according 
to the variable change operation associated to VS(.s'). Then, v 1 is normalized into a vector [v'], which is 
encoded into a word w G ±N{0, \} w . (In the case of multiple encodings, one of them can arbitrarily be 
chosen.) The word w corresponds to a single path of transitions leaving s, which is followed until a new 
implicit state s' is reached. Note that the states visited by this path between s and s' are explicit ones. 
The procedure then repeats itself from this state s' . 
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4.4.3 Examples 

An IRVA representing the set x\ > 1 A X2 < 2 A x\ — X2 < 1 in M 2 , considered in Figure [T^a), is given in 
Figure [2j Note that, since the set is not conical, the IRVA actually recognizes its representing cone, as 
discussed in Section 4.1 In this figure, implicit states are depicted by rounded boxes, and explicit ones 
by small circles. Doubled boxes represent in polarities. The vector spaces associated to implicit states are 
represented by one of their bases. Remark that the layout of the implicit states and the decision structures 
linking them closely matches the polyhedral components and their incidence relation as depicted in 
Figure [TJc), except for the initial state which corresponds to the apex of the representing cone. 

As an additional example, Figure [3] shows how the set x\ = 2 k X2 in M. 2 , discussed in the introduction 
of Section |3j is represented by an IRVA. In this case, the gain in conciseness is exponential with respect 
to RVA. 



5 Manipulation Algorithms 



5.1 Test of Membership 

A procedure for checking whether a given vector belongs to a polyhedron represented by an IRVA has 



already been outlined in Section 4.4.2 In the case of a polyhedron PCK" that is not conical, an 
IRVA can be obtained for its representing cone P C W +l , as discussed in Section 4.1 In this case, 
checking whether a vector (vi,V2, . . . ,v„) S W belongs to P simply reduces to determining whether 



(vi,V2, . . . ,v„, 1) belongs to P, which is done by the algorithm of Section 4.4.2 



5.2 Minimization 

An IRVA (n,Si,SE,so,A) can be minimized in order to reduce its number of implicit and explicit states. 
Since the transition relation A is acyclic, the explicit and implicit states can be processed in a bottom-up 
order, starting from the implicit states with the largest vector spaces. At each step, reduction rules are 
applied in order simplify the current structure. A first rule is aimed at merging states that are indis- 
tinguishable: If two explicit states share the same successors, they can be merged. In the case of two 
implicit states, one additionally has to check that their associated vector spaces are equal, and that their 
polarities match. The purpose of the second rule is to get rid of unnecessary decisions. Consider a state 
s (either implicit or explicit) with an outgoing transition that leads to an implicit state su representing a 
polyhedron component Q\. If all the implicit states Sj that are reachable from s are also reachable from 
si, then these implicit states represent polyhedral components Qj such that either <2; = Qi or Q\ -< Q\. 
The state s can then be absorbed into s\, provided that s is not an implicit state with a different polarity 
from the one of s\ . Note that this reduction rule correctly handles the case of a state s that is implicit and 
does not correspond to a polyhedral component, but to a proper subset of the component <2i represented 
by s\. For example, in M 2 , s may correspond to a unidimensional line x\ — X2 = covered by the larger 
universal component M? represented by s\. 

Property 1 Minimized IRVA are canonical up to isomorphism of their transition relation, and equality 
of the vector spaces associated to their implicit states. 

Proof sketch: The canonicity of a minimized IRVA srf representing a polyhedron P C W is the conse- 
quence of two properties. First, the minimization algorithm is able to identify and merge together implicit 
states that correspond to identical polyhedral components, as well as to remove the implicit states that 
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Figure 3: IRVA representing the set {(x\,X2) £ M 2 | x\ = 2 k X2\. 



do not represent such components. This yields a one-to-one relationship between the implicit states of 
si and the polyhedral components of P. Second, the transition structure leaving an explicit state s of s/ 



satisfies the following constraints. As discussed in Section 4.3 the state s corresponds to a component 
Q of P, and every prefix Wk of length k read from s defines a convex conical region R S Wl C W. If, in 
all sufficiently small neighborhoods of Q, the region R s ,w k covers a unique component Q' of P that is 
minimal with respect to the incidence order, then the path reading Wk from s leads to the implicit state s' 
corresponding to Q'. Provided that explicit states that have identical successors are merged, this property 
characterizes precisely the decision structure leaving s. Such structures will then be isomorphic in all 
minimized IRVA representing the same polyhedron. □ 

5.3 Boolean Combinations 

In order to apply a Boolean operator to two polyhedra Pi and P2 respectively represented by IRVA s/\ 
and s/2, one builds an IRVA s/ that simulates the concurrent behavior of &i\ and s^2- The procedure is 
analogous to the computation of the product of two finite-state automata. The initial implicit state of s/ 
is obtained by combining the initial states of s/\ and =2^2, which amounts to intersecting their associated 
vector spaces, and applying the appropriate Boolean operator to their polarities. Each time an implicit 
state s is added to s/, representing a polyhedron component Q, its successors are recursively explored. 



As explained in Section 4.3 each finite prefix Wk, of length k, read from s corresponds to a convex conical 
region R SWi C W. The idea is to check, in a sufficiently small neighborhood R of Q, whether R SjWk covers 
unique minimal components Q\ of Pi and Q2 of Pi, with respect to their respective incidence orders. In 
the positive case, one computes the intersection of the underlying vector spaces of Q\ and Q2. If the 
resulting vector space has a higher dimension than dim(VS(j)), as well as a non-empty intersection with 
R s .w k , a corresponding new implicit state is added to si '. In all other cases, the decision structure leaving 
s has to be further developed, which amounts to creating new explicit states and new transitions between 
them, in order to read prefixes longer than Wk- 

A key operation in the previous procedure is thus to compute, from an IRVA representing a polyhe- 
dron P, a component Q of P, and a given convex conical region C, the unique minimal component of P 
(if it exists) covered by C in the neighborhood of Q, with respect to the incidence order -<. This is done 
by exploring the IRVA starting from the implicit state representing Q. From a given implicit state s, the 
exploration only has to consider the paths labeled by words Wk such that C n Rs,w k 7^ ®> untn they reach 
another implicit state. Let S be the set of the implicit states reached this way. For each state in 5, one 
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checks whether its underlying vector space has a non-empty intersection with C. If this check succeeds 
for some nonempty subset of S, then the procedure returns its minimal component, or fails when such a 
component does not exist. Otherwise, it can be shown that the exploration can be continued from a sin- 
gle state chosen arbitrarily in S. The regions of space that are manipulated by this procedure are convex 
polyhedra, and can be handled by specific data structures Q. 

6 Conclusions 

We have introduced a data structure, the Implicit Real Vector Automaton (IRVA), that is expressive 
enough for representing arbitrary polyhedra in W, closed under Boolean operators, and reducible to a 
canonical form up to isomorphism. 

IRVA share some similarities with the data structure described in lfT4l . which also relies on decom- 
posing polyhedra into their components, and representing the incidence relation between them. The main 
original feature of our work is the decision structures that link each component to its incident ones, which 
are not limited to three spatial dimensions, and lead to a canonical representation. Furthermore, by imi- 
tating the behavior of RVA, we have managed to obtain a symbolic representation of polyhedra in which 
the membership of a vector can be decided by following a single automaton path, which is substantially 
more efficient that the procedure proposed in fPfl . 

The algorithms sketched in Section [5] are clearly polynomial. We have not yet precisely studied 
their worst-case complexity, since they depend on manipulations of convex polyhedra, the practical cost 
of which is expected to be significantly lower than their worst-case one. In order to assess the cost of 
building and handling IRVA in actual applications, a prototype implementation of those algorithms is 
under way. The example given in Figure [2]has been produced by this prototype. 

Future work will address other useful operations such as projection of polyhedra, conversions to 
and from other representations, and operations that are specific to symbolic state-space exploration al- 
gorithms. For this particular application, IRVA in their present form are still impractical, since they 
only provide efficient representations of polyhedra in spaces of small dimension. (Indeed, the size of 
an IRVA grows with the number of components of the polyhedron it represents, and simple polyhedra 
such as «-cubes have exponentially many components in the spatial dimension n.) We plan on tackling 
this problem by applying to IRVA the reduction techniques proposed in |4j, which seems feasible thanks 
to the acyclicity of their transition relation. This would improve substantially the efficiency of the data 
structure for large spatial dimensions. 
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